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Dihadron Fragmentation Functions (DiFF) encode information about the fragmentation pro¬ 
cess of a quark into a hadron pair (plus something else, undetected). Fragmentation has been hardly 
studied through models for pion pairs, the main knowledge about those DiFFs being for now fits 
from data. In particular, they can be extracted from the process [1,2], 

explored at Belle. 

We define the total P/, = Pi -|- P 2 and relative P = (Pi — P2)/2 momenta of the pair, with = 
Ml < G' = -<?" > 0 and q = k — k' the space-like momentum transferred and z is the sum of the 
fractional energies carried by the two final hadrons. The leading-twist cross section in collinear 
factorization, namely by integrating upon all transverse momenta but Rj- and Rj-, can be written 
as [3] 


- — - 1 -=-=- = ( 1+cos(0fi-h 0fi)Ae+e- ) , (1) 

d cos djdzdcos ddMijd^Rtizd cos OdMhd^a 4-n^ \ J 

where 62 is the angle in the lepton plane formed by the positron direction and Pj,, and the azimuthal 
angles and give the orientation of the planes containing the momenta of the pion pairs with 
respect to the lepton plane and the polar angle 6 which is the angle between the direction of the 
back-to-back emission in the center-of-mass (cm) frame of the two final hadrons, and the direction 
of P/i in the photon-proton cm frame. The da^ is the unpolarized cross section producing an 
azimuthally flat distribution of pion pairs coming from the fragmentation of unpolarized quarks, 
Di. In absence of data for the corresponding multiplicities, the unpolarized DiFFs have been 
parametrized to reproduce the two-pion yield of the PYTHIA event generator tuned to the Belle 
kinematics [1]. The term A^+g- represents the Artru-Collins asymmetry measured at Belle [4]. 
The latter is described in terms of chiral-odd fragmentation functions H^, which were therefore 
extracted from Belle data for the first time three years ago [1]. 

An update of the fit of has been recently proposed in Ref. [2] in which the error analysis 
is canied out using a Monte Carlo approach. This approach consists in creating A = 100 replicas 
of the data points. In each replica, the data point in the bin {zi,Mhj) is perturbated by a Gaussian 
noise with the same variance as the experimental measurement. Each replica, therefore, represents 
a possible outcome of an independent measurement. The fit has been performed with the choice of 
hadronic scale Go = IGeV^ and using two different values for as{Ml) = 0.125, 0.139 in the QCD 
evolution routine. In Fig. 1, we show the new results for the ratio 


R{z,Mh) 


|R| H<%z,M,r-,Ql) 
Mh D'{{z,Mh,Ql) ’ 


( 2 ) 


as a function of M/, for chosen values of z for as{Ml) = 0.139. The chisquare (per degree of 
freedom) distribution is centered around 1.6. The errorbands represent the 68 central replicas. 
There is only a mild dependence on the value of the strong coupling. No significative difference 
with the standad Hessian method of our previous extraction is noticed. 


1. Trans versify PDF 

Using (tt+tt^) SIDIS data off a transversely polarized proton target from HERMES [5] and 
proton and deuteron targets from COMPASS [6] together with the results for the DiFFs, a point- 
by-point extraction of transversity was performed for the first time in the collinear framework [7] 
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M,[GeV] 


Figure 1: The ratio R{z,Mh) as a function of M/, at = 1 GeV^ for three different z = 0.25 (shortest 
band), z = 0.45 (lower band at M/, ^1.2 GeV), and z = 0.65 (upper band at M/, ^1.2 GeV), with the value 
as{M 2 ) =0.139 used in the QCD evolution equations. 


and an independent parameterization of the valence components of up and down quarks proposed 
in Ref. [8]. An update of the extraction of the valence transversities, using data for identified pions 
from COMPASS [9], was released together with the new parameterization of the DiFFs [2]. The 
process of interest is Dihadron SIDIS off transversly polarized targets, 

l+N +H 1 +H 2 +X , (1.1) 


where I denotes the (unpolarized) lepton beam, N the nucleon target, H\ and H 2 the produced 
hadrons (here both are pions). The data range in x is restricted to x G [0.0065,0.29]. It allows one 
to access the (collinear) transversity, through the modulation from the azimuthal angle of the 
target polarization component St, transverse to both the virtual-photon and target momenta, and 
0R. The structure function attached to that modulation is written as the product of the PDF h\, and 
the DiFF [10], 


^UT 


{x,z,Mh\Q^) 


xY,e]h\{x\Q^) 


R| sinO 
Mh 




(1.2) 


Three different functional forms have been used for the parameterisation of the valence transver¬ 
sities, and their chisquare (per degree of freedom) distribution are centered in [1.4,2]. The main 
output, w.r.t the previous extraction, is the improvement of the acurracy of the fit in the mid-range 
X values for the up distribution, as reflected by the error band. The down quark valence distribution 
still tends to saturate the Soffer bound, a result that is somehow obtained in the single-pion SIDIS 
extraction with TMD evolution [11] as well, as is illustrated in Fig. 2. 

The tensor charge, defined as fhe firsf Mellin momenf of fhe valence fransversify PDF, is 
obfained by infegrafing over fhe physical supporf in x, 

SqviQ^) = tdxhfix,Q^) . (1.3) 

Jo 

Nolice lhal fhe error coming from fhe exlrapolalion of fhe PDF oulside fhe dala range is fhe main 
source of uncerlainly as shown in Table 3 of Ref. [2]. The value for fhe isovecfor fensor charge 
gT = 5uy — 5dy for fhe functional form relafed fo fhe so-called “flexible scenario" wifh as{M^) = 
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Figure 2: The up (left) and down (right) valence transversities as functions of x at = 2.4 GeV^. The 
darker band with solid borders in the foreground is our result in the flexible scenario with a^M^) = 0.125. 
The lighter band with dot-dashed borders in the background is the most recent transversity extraction from 
the Collins effect [12]. The central thick dashed line is the result of Ref. [11]. The thick blue lines indicate 
the Soffer bound. 

0.125 is gT = 0.81 ±0.44 at 2^ = 4 GeV^. It is in agreement with lattice determinations as well as 
with the other extractions from hadronic phenomenology [12, 13, 14], though the absolute value is 
slightly smaller than the lattice’s. 

2. Impact of Tensor Charge on BSM 

The determination of the tensor charge has an impact on beyond the standard model (BSM) 
physics searches through high precision beta decay observables. The latter allow us to probe cou¬ 
plings other than of the V —A type, which could appear at the low energy scale. Experiments using 
cold and ultra-cold neutrons, nuclei, and meson rare decays [15], are being performed, or have been 
planned, that can reach the per-mil level or even higher precision. These low energy measurements 
can be connected to BSM effects generated at TeV scales through effective field theories [16]. In 
this approach that complements collider searches, the new interactions are introduced in an effec¬ 
tive Lagrangian describing semi-leptonic transitions at the GeV scale including four-fermion terms, 
or operators up to dimension six for the scalar, tensor, pseudo-scalar, and V-i-A interactions, with 
effective coupling e,-, (/ = S, T,P,L,R) , e.g. [17]. 

The scalar (S) and tensor (T) operators contribute linearly to the beta decay parameters through 
their interference with the SM amplitude, and they are therefore more easily detectable. The ma¬ 
trix elements/transition amplitudes between neutron and proton states of all quark bilinear Lorentz 
structures in the effective Lagrangian which are relevant for beta decay observables, involve prod¬ 
ucts of the BSM couplings, e„ and the corTesponding hadronic charges, g,, i.e. considering only 
terms with left-handed neutrinos 


AXff=-Cspn-e{l-Y5)Ve-CT po^vH ■ - Y5)Ve , ( 2 . 1 ) 

where Cg = GfVudV^Esgs, and Cj = AGfVudy/lErgT^■ The tensor charge as determined here 
through dihadron SIDIS can be used to study the error on the BSM tensor effective coupling, 

*The latest global bounds on Ct are given in Refs. [18]. 
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Figure 3: Bounds on Cj- using various values of gj'- (turquoise-yellow) Lattice QCD; (green) Deeply virtual 
and rj production [13]; (blue) single pion SIDIS [12]; (red) dihadron SIDIS [2]. The dashed lines are 
for the future projections of the experimental extractions based on projections for JLab@ 12. The grey band 
gives the uncertainty in gj- assuming Agj- = 0 and the (weighted) average value for using the three lattice 
results plus the three hadronic results. 


currently still compatible with zero. In Ref. [19], the three hadronic extractions of the tensor charge 
— i.e. single-hadron SIDIS, dihadron SIDIS, 71° and rj in DVMP— have been considered to set an 
experimentally based limit on the possible new interaction, given the experimental bounds on Ct- 
Though the lattice simulations still give, in that particular case, the most strigent bounds, we have 
shown, see Fig. 3, that hadronic phenomenology will be competitive thanks to future experiments. 

As for future extractions, the dihadron SIDIS will be studied in CLAS12 at JLab on a proton 
target and in SoLID on a neutron target that will give together an improvement of ~ 10% in the 
ratio Agj/gT thanks to a wider kinematical coverage and better measurement of the d quarks 
contribution [20]. This estimation is based on a reanalysis of the fitting output with pseudo-data 
generated from projections with stastical and systematic errors for the mentionned experiments^. 
While the absolute value of the transversity can obviously not be determined yet, the magnitude of 
the predicted error is what matters here. 


3. Twist-3 PDF 

DiFFs also appear in asymmetries involving subleading objects. An extraction of the twist-3 
PDF, e{x), through the analysis of the preliminary data [21] for the sin^/;-moment of the beam- 
spin asymmetry for dihadron Semi-Inclusive DIS at CLAS at 6 GeV was proposed in Ref. [22]. 
Pion-pair production off unpolarized target in the DIS regime provide an access to the higher-twist 
Parton Distribution Functions e{x) and to Dihadron Fragmentation Functions. 

The e{x) PDF offers important insights into the physics of the largely-unexplored quark-gluon 
correlations, and its x-integral is related to the marginally-known scalar-charge of the nucleon, and 
to the pion-nucleon a-term. 

^See JLab proposals. 
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Figure 4; The extraction of the combination e^ = 4e“^ {xi, Qf) /9 — {xi,Qj) /9 in the WW scenario. The 

error bars correspond to the propagation of the experimental and DiFF errors. 

Beam Spin Asymmetry Alu presented in the proceedings [21], to leading-order in and 
leading term in the Partial Wave Analysis, is expressed in terms of FFs and collinear PDFs [23] 

ix,z,M,;Q\y) 

W{y) M |R| 4 

~~W)Q^r L,elf^{x-Q^)D\{z,M,-Q^) ’ 

(3.1) 

The dependence in {z,Mh) is factorized in the DiFFs and kinematical factors, leaving the depen¬ 
dence inX for the PDFs. The twist-2 functions are /i {x),Hf{z,Mh) andDi {z,Mh), while the twist-3 
functions are e{x) and G^{z,Mh). Under the assumption, related to the Wandzura-Wilczek approxi¬ 
mation, that the the twist-3 DiFFs are negligible, the PDF e{x) is directly extracted, point-by-point. 
The results are shown on Fig. 4. The average scale of the data is 2^ ~ 1.5 GeV^. A recent calcu¬ 
lation in the LFQM [24], in which only the kinematical part of the twist-3 PDF contributes, shows 
a good agreement with the present extraction, when evolved to 1.5 GeV^. The official publication 
of the data as well as complementary and future experiments will shed more light on the twist-3 
sector of the DiFF phenomenology. 

4. Conclusions 

We have presented various aspects of the rich phenomenology involving dihadron fragmenta¬ 
tion functions. Though dihadron SIDIS represent a (consequently smaller) subset of single-hadron 
SIDIS, the results achieved through these processes have been very promising so far. Future anal¬ 
yses — e.g. dihadron multiplicities at COMPASS— and experiments —especially in CLAS12 and 
Solid at JUab— will increase the data set, allowing for an improved knowledge on both DiFFs at 
leading and subleading-twist and collinear leading and subleading PDFs. 
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